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I.  INTRODUCTION 


Intermodal  dispersion  in  a  waveguide  results  in  widening  of  the  original  transmitted 
pulse.  Analogously,  multiple  arrivals  of  average  signal  energy  due  to  rough  surface 
scattering  result  in  dispersion  of  the  transmitted  pulse.  An  additional  source  of  pulse 
dispersion  in  a  waveguide  is  the  frequency  dependence  of  the  propagation  constant. 
Similarly,  frequency  dependence  of  the  propagation  constant  results  in  additional 
dispersion  of  the  time  dependent  scattered  field  --  especially  at  shallow  grazing  angles. 
By  assuming  an  arbitrary  narrowband  plane  wave  source  signal,  a  general  near-field 
expression  for  the  time  dependent  scattered  intensity  is  obtained  in  terms  of  statistics 
of  the  rough  surface  r-matrix,  and  equivalently,  in  terms  of  the  scattering  cross  section. 
Thus,  a  near-field  result  is  expressed  in  terms  of  a  far-field  entity,  the  scattering  cross 
section.  This  result  is  then  used  to  obtain  an  approximate  scattered  field  expression  for 
a  point  source.  These  results  apply  to  both  scattering  from  the  interface  as  well  as 
scattering  through  the  interface.  The  latter  case  is  of  interest  in  bottom  penetrating 
sonar. 

Although  the  theoretical  results  are  general  in  that  they  depend  only  on  the 
interface  7-matrix  statistics,  calculations  are  carried  out  using  first-order  perturbation 
theory  for  penetration  examples.  When  the  grazing  angle  of  the  incident  field  is  below 
the  critical  angle  in  relation  to  the  mean  surface,  only  the  zero-order  component  of  the 
transmitted  field  is  evanescent;  the  higher  order  components  contain  downward 
traveling  waves  (Moe  and  Jackson,  1994b).  Acoustic  intensity  as  a  function  of  depth  is 
computed  and  displayed  for  a  continuous  wave  (CW)  incident  plane  wave.  The 
resulting  profiles  of  intensity  versus  depth  are  compared  with  the  flat-surface  case.  The 
time-dependent  incoherent  intensity  is  computed  and  displayed  for  a  narrowband  point 
source.  These  computations  are  carried  out  using  a  low-frequency  cutoff  for  the  bottom 
relief  spectrum,  and  compared  to  the  flat-surface  case. 
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II.  SCATTERING  PROBLEM  BACKGROUND 


The  general  scattering  geometry  is  shown  in  Fig.  1.  An  arbitrary  pressure  field, 
\)/,(r),  is  incident  on  a  2-D  rough  interface  separating  the  homogenous  lossless  fluid  in 
medium  1  (water),  from  the  lossy  fluid  in  medium  2  (sediment).  The  zero-mean  rough 
surface  is  defined  by: 

z  =  hf{R),  (1) 

where  R  is  the  transverse  component  of  the  three-dimensional  position  vector  r ,  and  h 
is  the  root  mean  square  (RMS)  height  of  the  surface.  The  total  field,  \|ri(r),  above  the 
rough  surface  is  a  sum  of  an  arbitrary  incident  field,  \i/,(r),  and  the  resulting  scattered 
pressure  field,  v|/y^  (r): 

Vl(r)  =  Vi(r)  +  V/(r).  (2) 

The  boundary  conditions  consist  of  the  continuity  of  pressure: 

\iri=\K2.  (3) 

and  the  continuity  of  normal  velocity: 

—h  ■  V\|/i  =  —n  ■  V\)/2  (4) 

Pi  P2 

where  n  is  normal  to  the  surface,  pj  is  the  density  of  the  water,  and  P2  is  the  density  of 
the  sediment.  The  scattered  field,  v|/y  (r),  in  the  region  above  the  highest  point  on  the 

surface,  and  the  incident  field,  are  expressed  in  terms  of  the  following  (Weyl)  plane 
wave  expansions  (Devaney  and  Sherman,  1973): 

vj,^(r)  =  J  (5) 

and 

=  J  ,  (6) 
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where  K  is  the  two-dimensional  transverse  wave  vector  with  magnitude  K,  ki  \s  the 
wave  number  in  the  medium  above  the  surface,  and  itjpi,  is  the  z-component  of  the 
corresponding  three-dimensional  wave  vector,  with 

=  (7) 

For  the  continuous  wave  case,  a  factor  of  is  suppressed,  but  in  general,  the 
above  equations  are  the  Fourier  transforms  of  the  time  domain  field, 

¥l  (r)  =  Vi  (r,(D)  =  (8) 

The  field  below  the  rough  surface  is  expressed  as  a  plane  wave  expansion  of  waves 
traveling  in  the  negative  z  direction, 

xi/2_(r)  =  Jj2/i:^2-(K)e"'^P2(K)Zg/KR^ 


where, 

P2(K)  =  Vi-^^/A:| 


and 


h 


l  +  id 

V 


(10) 


(11) 


is  the  complex  wave  number  in  the  lossy  medium  below  the  surface  with  loss 
parameter,  6 .  The  speed  ratio,  v,  is  the  ratio  of  the  sound  speed  in  the  sediment  to  the 
sound  speed  in  the  water.  The  scattered  field  can  be  expressed  in  terms  of  the 
incident  field  and  the  T-matrix,  Tif 

^f{K)  =  jd^K'Tu{K,K')'i'i{K')  (12) 

Likewise,  the  downward  component  of  the  scattered  field  penetrating  below  the  surface 
is  expressed  in  terms  of  the  T-  matrix,  rj2: 
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4'2_(K)  =  Jrf^*:'7-,2(K,K')4',.(K').  (13) 

The  scattered  field  for  the  case  of  an  incident  plane  wave  is  found  by  substituting, 

>i'i(K)  =  5(K-K,.),  (14) 

into  Eq.  (12),  and  the  result  into  Eq.  (5),  yielding 

\\ff  (r)  =  J (15) 

The  scattered  field  penetrating  the  surface  for  the  incident  plane  wave  case  is  found  in 
the  same  way  from  Eqs.  (9),  (13)  and  (14), 

X|/2_(r)  =  J^/2ii:ri2(K,K,)e“'*^2p2(K)ZgfKR  (^5) 


III.  GENERAL  RESULT  FOR  INCIDENT  PLANE  WAVE 

The  purpose  of  this  section  is  to  derive  an  expression  for  the  average  intensity  time 
series  at  the  receiver  position,  r,  due  to  a  narrowband  plane  wave  source.  Although 
the  following  derivation  considers  the  more  general  case  of  scattering  through  the 
rough  interface,  it  also  applies  to  scattering  from  a  rough  interface. 

A.  Time  domain  received  signal 

A  narrowband  source  signal  s{t),  with  Fourier  transform: 

oo 

5(o))=  J^r5(r)e'“',  (17) 

— oo 

can  be  represented  as: 

s(<)  =  Re{«(()e-'®^},  (18) 
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where  u{t)  is  the  corresponding  complex  baseband  input  signal,  with  Fourier  transform, 
C/((d),  and  is  the  carrier  frequency,  or  center  frequency  of  the  pulse.  An  incident 
plane  wave  pulse  with  direction  denoted  by  the  transverse  unit  vector,  a,-,  and  incident 
grazing  angle  9,  ,  is  expressed  in  terms  of  its  2-D  Fourier  transform  as: 

=  ©>0,  (19) 

where 


K/ =— cos0,a/,  (20) 

is  the  incident  transverse  wave  vector,  and  Zs  is  taken  to  be  an  arbitrary  source  height. 
Using  Eqs.  (9),  (13),  (17),  (18),  (19),  and  chapter  3,  Proakis  (1989),  the  baseband 
pressure  signal  (r,r)  at  a  position  r  =  (R,z)  below  the  surface  is  expressed  as: 


(r,0  =  j  d^K  r,2(K,K,.,®'+a,,) 


^^jsine,z^((o'+co,:)/ci^-rKP2(i^,co'+(Oc)z(co'+(Oc)/qgjKR 


where 

k  =  (22) 

is  the  complex  wave  number  ratio, 

K  ■  =  ^  cos  6,a/ ,  (23) 

Cl 

and  we  have  explicitly  shown  the  dependence  of  p2(K)  on  frequency  as  well  as  on  the 
magnitude  of  K. 


B.  Time  dependent  intensity 

The  time  dependent  intensity  at  position  r  can  be  expressed  as: 
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(24) 


Combining  Eqs.  (21)  and  (24)  yields: 

/2(r,/)  =  J  d<Si'\  d(S)"\u{(ii')U\ai")  e"'®''  e*®"' 

xj  d\"]d\'  ^7i2(K',K<,co'  +  co  Jri2(K",K<',o)"+a)^)) 
^  giK'Rg-iK"-RgJ  sin  Gj-z^co'/ci  g-i  sin  e,z^(o7ci 

^^-iKp2(A:',(o'+CDc)zco7c,^i(Kp2(^:",co"+(Oc))*z(a)"+cOc)/ql 


(25) 


where  the  transverse  wave  vectors  for  the  incident  field  at  the  angular  frequencies  co", 
and  CD'  are 


cos6,a,  ,  =^^-^^cos9,d,-. 

q  q 


(26) 


Define  Ci2(K',K",K;,Kf,cOc+co',(Oc  +  co"),  (see  for  example  Zipfel  and  DeSanto,  1972; 
Voronovich,  1995)  such  that 


(ri2(K',K;,(0'+C0,)7i2(K",Kf,(0"  +  0),)) 

=  Q  2  (K',  K",  k;  ,  K;',  0c  +  0)',  0c  +  0")5(K"  -K'  +  K'i-K'f) 


(27) 


In  the  above  expression,  the  subscript  is  used  to  represent  scattering  into  the  sediment 
and  the  dependence  of  T12  on  0  is  included  in  the  argument.  Substituting  Eq.  (27)  and 


Kf-K<  =  - — ^cose^a,-  =  Kj 
Cl 

into  Eq.  (25),  and  changing  the  integration  variable  K'  to 


K  =  K'  +  K^/2, 


(28) 


(29) 


yields: 
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(30) 


2 

/2(r,0  =  J  Jco'j Jco"l/(a)')t/*(a)") 

xJj%Q2(K-K^/2,K+Kj/2,K;,K;',o)^  +  {o',a)c  +  co")e“*^-*^ 

^  gi  si  n  QjZs  co'/q  g-i  sin  0/z^  ©"/c) 

^^-/Kp2(lK-Kj/2|,©c+©')z©7q  /(K|i2(|K+Krf/2|,©,+©''))*z(©c+©")/q 
X,  c  c 

Equation  (30)  is  now  evaluateid  by  making  a  few  approximations.  For  the  case  of  a 
narrowband  source  signal,  we  can  assume  Q2  '®  ®  slowly  varying  function  of 
frequency, 

Q2(K-K^/2,K  +  Krf/2,K;,Kf,(l)^  +  0)",0)^  +  C0')  =  q2(K.K,Ki,Kj,<0<,,0)J 

sq2(K,K,.m,) 

where  the  relation  (Thorsos  and  Jackson,  1989;  Berman,  1992), 

Q2(K,K,-,(o^)5(K-K'0  =  (ri2(K,K,-,coJri2(K^KpCO,))  (32) 

is  a  special  case  of  Eq.  (27),  and  K,  is  evaluated  at  the  center  frequency: 

K,- =-^cose,d,.  (33) 

First-order  perturbation  theory  can  be  used  to  show  that  Equation  (31)  is  a  valid 
approximation.  Also,  define  the  vertical  component  of  the  wave  vector  in  medium  2, 

B2{K,(i)c+Oi',(Oc+  co")  s  ((o)^  +  (d')/ci  )kP2  (jK  -  K  j  /2|,  ©c  •+©'),  (34) 

Analytical  expressions  for  dispersion  in  a  waveguide  can  be  found  by  expanding  the 
exponent  of  the  time  dependent  field  in  a  power  series  (for  example,  see  Ishimaru, 
1991).  Since  Kj  =  \Kj\  is  a  function  of  both  ©'  and  ©",  (©7ci)kP2  is  expanded  in  a 
power  series  in  ©'  and  ©",  for  the  case  of  a  narrowband  signal; 
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(35) 


^2  co^  +  a)',(Oc  +  (o")  = +1^ 


1 

+— 

2 


3  B2  ,2  ^  ^  ^2 

- ^0)'^  +  2— — ^ 

0(o'^  00)9© 


©  ©  +-  ^ 


0©' 


where  the  derivatives  are  evaluated  at  ©'  =  ©"  =  0,  and 


(36) 


where  the  wave  number,  fcj,  in  the  above  equation,  and  all  following  equations,  is 
evaluated  at  the  center  frequency  (^j  =  ©^/q ).  Solving  for  the  derivatives  in  Eq.  (35) 
and  combining  like  powers  of  ©'  and  ©"  in  Eq.  (30)  yields, 

/2(r,r)  =  J  d^K  |Q2(K,K,-,©^)e^^i^^'"^''^2(^.®c))2 

x_L  r^(o"t/*(©")e'“"('-'i-'2re'“"^^V4  ,  (37) 

271'’ 

X  ±  j  da'U{a')  -'2)e-“'“"2r''“''  '’/“j 


where 


h 


Rg; 


cosSj 


+— sin6/  +  z 


&;■  •Kcos6,- 
q^i 


r 

Re 

V 


Kp2(/^,C0,)J 


and 


-ZK/q 

^  P2(^2“c)’ 


(38) 


(39) 


The  physical  significance  of  these  parameters  will  be  discussed  later.  Note  that  ti  is 
real  and  t2  is  complex.  The  complex  coefficient 

P  =  {(2k^/ )(«/•  K  COc  cos  e,  /  q^  -  (/:/q  )^  )  -  z  Im(l/  B2c  )(cos  6^  /q  )^ 

\S  .  (40) 

-zlm^l/  P2c)(®i  •Kcos9j/q)^|z 
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and  the  real  coefficient 


Q  =  Im|(d/  •  K  cos  9,-  / B^c  )(k^C0c / c?  -  dj  ■  K  cos  9,-  j {icf ))  -  cos^  9,/ (2cf  B2c )}^  (41 ) 

lead  to  pulse  dispersion.  Defining 
j2(,)  =  J-  f 

.  (42) 

271  •’ 

results  in  a  simple  expression  for  the  time  dependent  intensity  in  medium  2: 

In  addition  to  being  a  function  of  the  complex  argument,  t ,  g{t)  is  a  function  of  K,  the 
incident  field  direction,  and  fluid-sediment  parameters.  Note  that 

U=h'^h-  (44) 

is  also  a  complex  function  of  K.  For  an  arbitrarily  small  pulse  bandwidth  (CW),  the 
dispersion  terms  are  arbitrarily  small,  and  g{t)^  u{t). 

While  Eq.  (43)  involves  an  integral  over  the  transverse  wave  vector,  it  can  be 
converted  to  an  integral  over  the  scattering  surface  through  the  change  of  variable, 

R'  =  R-vr^K//ti.  (45) 

As  shown  in  Fig.  2,  R'  is  the  location  of  a  small  scattering  patch  on  the  surface,  and 
rj  =  |rj|  is  the  distance  from  this  patch  to  the  field  point,  r , 

rj  =  >-Rf +7^  (46) 

As  the  variables  R'^  and  /?'  range  over  all  real  values,  Eq.  (45)  constrains  K  to  the 
range  K<ki/v.  This  change  of  variables  forces  neglect  of  evanescent  waves. 

The  Jacobian  follows  from  Eq.  (45), 
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(47) 


Jd+z^  xdyd 

Xdyd  Xd+z^ 

-  -2 
“  2  4^’ 

where  x^,  and  yj  are  the  transverse  coordinates  of  r^.  The  geometric  significance  of 
these  definitions  is  apparent  when  one  assumes  that  the  loss  parameter  is  small 
(5«l),  Then 


^lKP2{/i:,03^)  =  VKf-^ 

=  (^j  sin  02  /  v)^l  +  ib/  sin^  02  j 


(48) 


where 

sin  02  =  =  \z\/rj  (49) 

is  the  sine  of  the  scattered  grazing  angle  defined  in  Fig.  2.  Dropping  terms  of  second 
order  and  higher  in  the  loss  parameter,  the  significance  of  rj  and  t2  becomes  apparent, 
as 

?!  =  R'  d,' cos0j7ci sin0,/ci  (50) 

is  the  time  required  for  the  incident  plane  wave  front  with  direction  (d,  cos0,-,sin0,)  to 
travel  from  the  source  reference  point,  (0,0,2^)  to  the  scattering  patch,  (R',0),  and 

^^{h)  =  fd/^2  (^‘•) 

is  the  time  required  for  a  spherical  wave  scattered  from  the  patch  at  (R',0)  to  travel 
through  the  sediment  to  the  field  point  (R,z),  z<0.  According  to  Eq.  (48),  the 
absorption  exponent  in  medium  2  is  given  by 

1  c 

Im(t,K|32(K,(0^))lz|  =  -^|z|/sin  02  =  Jm{k2)rj. 


(52) 


Substituting  Eqs.  (47)  and  (52)  into  Eq.  (43)  results  in: 


'  '  V  rj 

Substituting  the  following  quantity, 

Oi2(a,a,)  =  M  q^(k,k,,cOc), 

V 

a  =  KIK,  and  a,  =K,  //i:,- 


(53) 


(54) 

(55) 


into  Eq.  (53)  yields, 

l2(r,t)  =  \d^R'- 


(56) 


The  interpretation  of  Ci2  as  a  scattering  cross  section  will  be  justified  later  in  this 
section.  Since  K  =  /:](R-R')/(vrj),  t2,  g,  and  therefore  P,  and  Q  are  functions  of 

|R-R'l.  Limiting  the  integration  in  Eq.  (43)  to  the  non-evanescent  regions,  K<kifv, 
noticeably  affects  the  integral  only  when  within  a  wavelength  of  the  surface.  Note  that 
the  CW  result  follows  from  Eq.  (56)  by  simply  setting  g{t-td)  =  ^-  For  this  case  of  an 
arbitrarily  narrowband  source  signal,  the  result  in  Eq.  (43)  is  exact,  and  Eq.  (56)  is 
exact  everywhere,  except  very  close  to  the  surface. 

Although  the  intensity  of  the  field  scattered  into  the  water,  /i(r,r),  is  similarly 
derived,  the  corresponding  expression  for  /i(r,r)  is  most  easily  obtained  as  a  special 
case  of  Eq.  (56)  by  substituting  medium  1  parameters  for  medium  2  parameters.  For 
example,  ki  is  substituted  for  k2,  Pi  for  KP2,  etc.  For  this  special  case,  Eq.  (56) 
simplifies  to 


/i(r,r) 


1“ 


^R' 


8h{^-td)< 


(57) 
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where  from  Eq.  (54) 


an(a,a,)  =  k}  sin^(e2)|Cii(K,K,-,a)c)p.  (58) 

In  this  case,  62  is  also  the  scattered  field  grazing  angle  measured  from  the  mean 
horizontal  plane,  |zl/r^  =sin(02).  Again,  the  quantity  Qj  is  found  from  the  relation  (Eq. 
32) 

Ql(K,K^(O,)8(K-K'0^{7ii(K,K,-,(D,)ri*i(K'',K,-,co,))  (59) 


There  is  no  loss  in  this  case  --<2  =  0,  and  therefore, 

—  f  rfco'  f/((o')e"'“''e"'®'^  ^ 

2n^ 


sn(')“ 


P]  1  E  cos9,*  cos»t'/ cf  -(AT/c,)^  J, 


(60) 

(61) 


and 


(62) 


Note  that  is  still  given  by  Eq.  (44),  but  the  propagation  time  from  the  surface  patch  to 
the  receiver,  or  observation  point  above  the  surface  is  a  real  quantity  — 


t 

2  pi(i^,co,)' 


(63) 


As  shown  in  Voronovich  (1995),  the  quantity  Oii(a,d,)  in  Eq.  (58)  is  actually  the 
scattering  cross  section  per  unit  area  defined  in  Ishimaru,  (1978). 


Oii(d,d,)  = 


r-r 


a|2 


M 


(64) 


where  v)/,  is  a  plane-wave  field  incident  on  a  surface  patch  of  area  AA,  and  d,-  is  the 
unit  vector  in  the  direction  of  propagation  of  the  incident  field;  \|iy  is  the  scattered  field 
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at  a  long  distance  |r-r'|  from  the  surface,  In  the  direction  denoted  by  the  unit  vector, 
a.  Similarly,  Ci2,  is  the  scattering  cross  section  per  unit  area  of  a  rough  surface 
relating  the  incident  field  intensity  to  the  scattered  field  intensity  in  medium  2  (the 
sediment),  defined  as: 


ai2(d,a,)  = 


^2Iin(/:2)|r-r'| 


(65) 


where  is  the  field  at  the  receiver  in  the  sediment,  a  distance  jr-r'l  from  the 
surface  patch,  and  the  unit  vector  a  represents  the  transverse  direction  of  the  scattered 
field  in  the  lower  medium.  The  result  given  in  Eq.  (56)  (or  Eq.  57)  is  remarkable,  as  the 
cross  section,  Eq.  (65),  is  a  far-field  entity.  Close  inspection  of  Eq.  (56)  shows  that  it  is 
the  "naive"  sonar  equation,  when  u{t)  is  substituted  for  g(r),  in  which  the  scattering 
intensity  is  obtained  by  integrating  the  bistatic  cross  section  over  the  interface  with 
appropriate  attenuation  due  to  spreading  and  sediment  loss.  Equations  (32)  and  (54) 
are  a  convenient  way  to  find  Cu  in  theoretical  developments;  Eqs.  (58)  and  (59)  are  a 
convenient  way  to  find  Oj  i . 


C.  Dispersion  of  a  Gaussian  input  pulse 

Two  types  of  pulse  dispersion  occur  in  this  problem.  One  type  is  simply  due  to 
scattering  and  is  treated  by  the  integral  in  Eq.  (56)  over  the  scattering  surface.  The 
other  type  of  dispersion  is  due  to  the  frequency  dependence  of  the  propagation 
constant,  and  is  embodied  in  Eqs.  (42).  This  subsection  focuses  on  this  form  of 
dispersion.  Consider  the  case  of  an  input  pulse  with  a  Gaussian  envelope, 

u{t)  =  ,  (66) 

with  Fourier  transform 
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(67) 


Substituting  Eq.  (67)  into  Eq.  (42)  and  integrating  over  co'  and  co"  yields: 
2,.  tj  -t'^lq'2  -(t*+2tQlq'^f  jq"^ 

q  q 

where 

q'^yjtj  +  iP , 


(68) 


(69) 

(70) 


and 


P"  =  P*-i{2Qjqf, 


(71) 


Substituting  Eqs.  (69)  and  (70)  into  Eq.  (68),  and  simplifying,  yields  an  expression 
equivalent  to  Eq.  (68) 


(72) 


As  mentioned,  Re(r^)  =  Re(/2)+ti  represents  the  propagation  delay  of  the  pulse.  The 
imaginary  part  of  t2  is  a  result  of  modeling  the  sediment  as  lossy.  Since  higher 
frequency  signals  are  attenuated  more  than  lower  frequency  signals,  the  lower 
frequency  components  of  a  narrowband  signal  will  be  less  attenuated  than  the  center 
frequency.  The  imaginary  component  of  t2,  along  with  Q,  and  the  imaginary 
component  of  P  compensate  for  typically  excessive  loss  in  the  CW  loss  term  —  the 
exponential  term  in  Eq.  (43)  or  Eq.  (56).  In  addition,  P  is  a  frequency  dispersion  term 
that  results  in  pulse  broadening  of  the  signal  below  the  rough  surface  (specifically  the 
real  component  of  P).  Note  that  for  the  zero  loss  case,  the  parameter  <2  =  0.  The 
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effect  of  loss  and  the  parameters  P  and  Q  is  seen  in  Fig.  3  for  g{t-td)  from  Eq  (72) 
evaluated  at 


=  sH-‘2i)'  '2i  =  Im(<d)  =  Im('2) 


l=Re(lj) 


(73) 


and  plotted  as  a  function  of  scattered  grazing  angle  for  lossy  and  zero-loss  sediment. 
In  Fig.  3,  the  pulse  length  parameter,  is  equal  to  two  periods  of  the  center,  or  carrier 

frequency  for  both  cases  shown.  When  the  loss  and  dispersion  are  small,  |g(-r2/)| 
approaches  m(0)  =  1.  The  smaller  the  scattered  grazing  angle,  the  greater  the 
dispersion,  and  the  smaller  one  would  expect  the  pulse  peak  amplitude.  However,  the 

loss  of  the  CW  signal  represented  by  the  term,  1^  jg  loss  of 

the  narrowband  signal,  and  the  propagation  distance  becomes  large  for  small  scattered 
grazing  angles  (z  is  fixed).  Even  though  there  is  dispersion,  the  loss  for  the  narrow 
band  case  is  significantly  less  than  the  CW  case,  and  |g(-/2,)|>l  •  In  the  case  of  Fig. 

3,  with  depth  of  0.3  m,  |g(-r2()|  is  greater  than  one  for  very  small  grazing  angles  since 
the  propagation  distance  becomes  very  large.  An  example  with  zero  loss  always  has 
P"  =  P* ,  t2i  =  0 ,  and  it  is  always  true  that  |g(0)|  <  1 . 


IV.  ANALYTICAL  SCATTERING  MODEL  (PERTURBATION  THEORY) 

In  the  examples  of  the  following  section,  the  intensity  of  the  field  in  medium  2  is 
divided  into  coherent  and  incoherent  parts.  The  incoherent  intensity  is  treated  using 
first-order  perturbation  theory,  and  the  coherent  intensity  is  treated  in  zeroth  order  -  the 
flat  surface  solution  is  used.  There  is  a  slight  inconsistency  in  this  approach,  in  that  the 
incoherent  intensity  is  proportional  to  the  second  power  of  k^h,  while  the  coherent 
intensity  is  only  computed  to  an  accuracy  of  zeroth  order.  As  Rice  (1951)  has  shown. 


16 
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the  coherent  intensity  to  second  order  in  k^h  is  obtained  by  subtracting  the  power 

carried  by  the  incoherent  intensity.  For  small  k^h,  this  correction  is  necessarily  slight. 

A.  Incoherent  Intensity  using  the  first-order  perturbation  approximation 

In  this  section,  the  first-order  perturbation  results  are  given  for  the  scattering  cross 
section  for  insertion  in  Eq.  (56).  The  first-order  r-matrix  is  proportional  to  the  2-D 

Fourier  transform  of  the  surface,  F(K),  (see  the  Appendix)  and  is  given  by: 

7i2(K,Ki)  =  (*,;.)rW(K,K, .)  =  (»:, *)Hi2(K,K,)F(K-K,), 

(74) 

where 

(75) 

is  a  function  of  the  flat-surface  reflection  coefficient,  r,  the  flat-surface  transmission 

coefficient. 

7'(K)=i-Hr(K) 

(76) 

as  well  as  the  factors  a(K",K)  and  b(K'\K)  defined  in  the  appendix.  Using 

h^{F{K"  -  K,)F(K'  -  K,' )}  =  W(K'  -  K,)5(K"  -  K') 

(77) 

together  with  Eq.  (32),  and  Eq.  (74)  yields: 

q2(K,Ki,o,,)  =  |t,Hi2(K,Ki)fw'(K-K,-). 

(78) 

An  expression  for  the  first-order  incoherent  intensity  results  from  substituting  Eq, 
into  Eq.  (43).  Equivalently,  the  first-order  intensity  is  found  from  Eq.  (56),  where 

(78) 

og  (dy- ,  a, )  =  .K,  1  V(Ky- -  K, ) , 

(79) 
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the  scattering  cross  section  from  first-order  perturbation  theory,  is  found  from 
substituting  Eq.  (79)  into  Eq.  (54).  Although  the  cross  section  in  this  approximation  is 
actually  second  order  in  k^h,  (Thorsos  and  Jackson,  1989),  we  will  refer  to  it  as  the 
first-order  cross  section,  and  will  refer  to  the  corresponding  scattered  intensity  as  the 
first-order  intensity. 

B.  Zero-order  intensity 

The  zero-order,  or  flat-surface  transmitted  coherent  intensity  for  an  incident  plane 
wave,  is  given  by  the  magnitude  squared  of  Eq.  (16),  with  r/2^(K,K,)  given  by  (A6): 


C.  Roughness  spectrum 

The  two  dimensional  Gaussian  random  process  /(R)  describing  the  seafloor 
surface  is  assumed  to  be  isotropic,  with  a  roughness  spectrum,  W(K),  in  the  form  of  a 
filtered  power-law  (Moe  and  Jackson,  1994b), 

=  ,81) 

with  RMS  height. 


Here,  T  refers  to  the  gamma  function  --  not  the  reflection  coefficient.  The  factor 

is  obtained  by  subtracting  a  Gaussian  weighted  moving  average 

y 

(Papoulis,  1984)  of  the  surface  from  itself. 
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V.  APPLICATIONS  AND  CALCULATIONS 


Although  the  theoretical  results  in  this  paper  are  general  and  include  scattering 
back  into  the  water,  the  examples  deal  with  scattering  through  a  rough  fluid-fluid 
interface  (z  <  0).  Both  plane-wave  and  point  sources  are  considered,  and  the  scattered 
intensity  is  computed  for  the  flat-surface  case. 

A.  Continuous  plane  wave  source 

The  first-order  intensity  of  the  field  penetrating  through  a  rough  surface,  with 
roughness  spectrum  given  in  Eq.  (82),  due  to  a  CW  incident  plane  wave  is  calculated 
using  Eqs.  (78)  in  Eq.  (43),  and  setting  =  Using  parameters  appropriate  to  a 

sandy  seafloor,  the  incoherent  acoustic  intensity  is  computed  as  a  function  of  depth. 
Since  the  grazing  angle  in  this  example  is  below  critical,  the  zero-order  (flat-surface) 
field  is  evanescent.  Near  the  surface,  the  incoherent  intensity  decays  with  increasing 
depth  at  rate  comparable  to  the  zero-order  coherent  intensity,  but  decays  at  a  rate 
dictated  by  the  sediment  loss  further  from  the  surface.  This  example  (Fig.  4)  shows  that 
roughness  has  a  significant  effect  on  the  field  transmitted  through  a  surface. 


B.  Incoherent  field  due  to  a  point  source  including  dispersion 


Expressing  the  2-D  integral  in  Eq.  (56)  in  terms  of  a  summation  over  the  surface 
yields 


/jCr.OsAAj 


rin 


(83) 


where 

r2n=\r-K 


(84) 
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is  the  distance  from  surface  patch  n ,  at  position  r;^  of  size  AA  to  the  observed  field 
point,  r ,  and 


R-R' 

R-R' 


(85) 


The  propagation  time  from  this  surface  patch  to  the  receiver,  Re(r2n).  is  also  given  by: 


R«('2/i)  =  ''2n/<^2' 


(86) 


As  the  area  of  each  surface  patch,  AA,  approaches  zero,  the  summation  in  Eq.  (83) 
approaches  the  exact  integral  in  Eq.  (56). 

An  expression  for  the  incoherent  first-order  field  below  the  surface  due  to  a  point 
source  follows  from  Eq.  (83).  Here,  the  surface  patch  (area  AA)  is  chosen  sufficiently 
small  in  relation  to  the  distance  from  the  point  source  to  surface  patch  n ,  to  insure 
that  the  field  incident  on  the  surface  is  a  plane  wave.  The  scattered  or  penetrating  field 
intensity  is  given  by: 


l2(r,t)  =  AA^ 


2  2 
nnrin 


-g 


(87) 


where, 


K-^s 


(88) 


and  Rj  is  the  transverse  coordinate  vector  of  the  source.  The  total  propagation  time  is 
given  by  the  real  part  of: 


(89) 


with 


(90) 
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When  the  approximation  g{t)-u{t)  is  made,  Eq.  (87)  is  in  the  form  of  a  convolution, 
reducing  the  computation  time  significantly. 

Using  the  expression  for  the  first-order  scattering  cross-section  given  in  Eq.  (79), 
the  approximate  incoherent  intensity  pulse  in  lossy  sediment  below  a  rough  surface 
due  to  a  Gaussian  narrowband  point  source  above  the  surface  is  found  for  two  cases. 
In  Fig.  5a,  the  incoherent  pulse  time  series  at  a  position  below  a  rough  surface  in  a 
lossy  medium  due  to  a  point  source  is  found  using  Eq.  (72)  for  g^(r)  in  Eq.  (87),  and  is 
compared  to  the  result  obtained  setting  the  dispersion  terms  P  and  Q  to  zero  — 
g{t)  =  u{t).  Although  the  resulting  pulse  shapes  are  close  in  magnitude,  as  well  as  peak 
arrival  time,  including  the  dispersion  terms  does  result  in  an  incoherent  intensity  pulse 
that  has  a  smaller  peak  magnitude,  and  is  wider  than  the  approximation  obtained  by 
setting  g{t)  =  u{t).  In  this  example,  the  parameter  in  Eq.  (66)  is  equal  to  two  cycles  of 
the  center  frequency  of  =20Khz.  This  frequency  results  in  kih  =  0.25,  which  is  within 
the  region  of  accuracy  for  first-order  perturbation  theory  (Thorsos  and  Jackson,  1989; 
Thorsos,  1990).  In  Fig.  5b,  the  point  source  is  higher,  but  the  incident  angle  is  the 
same.  In  this  example,  including  the  dispersion  parameters  has  a  slightly  greater 
effect. 


C.  Zero-order  calculations 


Setting 

'F2-(K)  =  9-(K)'F,(K) 

and 

\i/  _ _ i _ 

2jrt,P,(K) 


(91) 


(92) 


in  Eq.  (9)  results  in  an  expression  for  the  field  penetrating  a  flat  fluid-fluid  interface  due 
to  a  point  source  at  (Ri,Zi)  with  unit  magnitude  at  Im  from  the  source.  The  resulting 
integral  simplifies  to: 

\}/2-(r)  =  /  j  dK 

0  «lPlW 

where  Jq  is  the  Bessel's  function  of  order  zero.  For  comparison  purposes,  the  zero- 
order,  or  flat-surface  coherent  intensity  is  plotted  in  Fig.  6a  along  with  the  first-order 
intensity  pulse  using  the  same  parameters  as  in  Fig  5a.  In  this  example,  the  incident 
field  is  below  the  critical  angle.  Although  the  zero-order  component  can  be 
approximated  from  Eq.  (93)  using  the  method  of  stationary  phase  to  solve  the  resulting 
integral  (see  for  example  Brekhovskikh,  1980;  Ishimaru,  1991)  when  the  incident  field  is 
above  the  critical  angle,  a  numerical  method  is  more  suitable  when  the  incident  grazing 
angle  is  close  to  critical  (Westwood,  1989).  The  incident  grazing  angle  is  defined  here 
to  be  the  angle  between  the  mean  surface  and  the  line  containing  the  source  and 
receiver  points.  Figure  6a  shows  an  “exact”  solution  of  the  flat-surface  intensity  time 
series  pulse,  along  with  the  first-order  pulse  (including  dispersion  terms)  from  Fig.  5a. 
For  this  example,  the  zero-order  intensity  pulse  arrives  earlier  and  with  greater 
magnitude  than  the  first-order  pulse.  In  Fig.  6b,  the  source  is  moved  further  from  the 
surface  than  the  previous  two  examples.  Even  though  the  sediment  parameters, 
receiver  point,  and  incident  angle  are  unchanged  from  Fig.  6a,  the  incoherent  first- 
order  intensity  pulse  is  greater  in  magnitude  than  the  zero-order  intensity  pulse.  As  the 
source  is  moved  further  from  the  surface,  the  stationary  phase  path  for  the  zero  order 
pulse  is  longer,  and  the  zero-order  pulse  amplitude  is  more  attenuated  in  the  lossy 
sediment.  Since  Fig.  4  is  a  plane  wave  CW  version  of  this  same  example,  one  would 
expect  that  the  first-order  intensity  would  be  greater  than  the  zero-order  intensity  for 
the  source  point  sufficiently  far  from  the  surface. 
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VI.  CONCLUSION 


A  general  analytical  expression  for  the  time  dependent  narrowband  scattered 
intensity  from  and  through  a  fluid-fluid  rough  interface  can  be  expressed  in  terms  of  the 
second  moment  of  the  r-matrix,  or  equivalently,  in  terms  of  the  scattering  cross  section 
per  unit  area.  For  the  special  case  of  a  CW  incident  plane  wave,  this  result  is  exact, 
which  is  suprising  as  the  scattering  cross  section  is  a  far-field  entity.  Using  first-order 
perturbation  theory  to  obtain  the  scattering  cross  section  per  unit  area  yields  a  simple 
expression  for  the  incoherent  scattered  time  dependent  intensity  pulse.  This  expression 
is  especially  tractable  when  the  source  pulse  is  Gaussian. 

The  above  formalism  applied  to  the  scattering  through  a  rough  surface  shows  the 
effect  of  roughness  on  penetration  through  a  surface.  When  the  grazing  angle  of  a 
plane  wave  incident  field  is  below  the  critical  angle  in  relation  to  the  mean  surface,  the 
zero-order  component  of  the  transmitted  field  is  evanescent,  and  does  not  penetrate 
deeply  into  the  seafloor.  Higher  order  components  contain  downward  traveling  waves, 
which  can  increase  the  depth  of  penetration  of  sound  relative  to  the  flat-surface  case. 
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APPENDIX:  PERTURBATION  THEORY  RESULTS 


The  fields  and  r-matrices  are  expanded  in  their  zero-order  and  first-order 
components,  with  the  factor  k^h  displayed  explicitly; 


'F(K)  =  fci 

(A1) 

T(K,Ki)  =  T^^\K,Ki)  +  kihT^^\K,Ki) 

(A2) 

Zero-order 

The  zero-order  scattered  field,  or  reflected  field: 

¥f\K)  =  r(K)'P,(K) 

(A3) 

where  r(K)  is  the  flat-surface  reflection  coefficient. 

,  P|3,(K)-k|32{K) 

^  P|3i(K)  +  k|32(K) 

(A4) 

Comparing  Eq.  (12)  to  (A3),  the  reflection  coefficient 

matrix  by: 

is  related  to  the  zero-order  T- 

r,<'»(K,Ki)  =  r(K)5(K-K,) 

(A5). 

The  transmitted  field  is  given  in  terms  of  the  incident  field  and  transmission  coefficient, 

'F^°^(K)  =  T(K)'P,(K)  (A6) 

Comparing  Eq.  (13)  with  (A6),  the  transmission  coefficient  is  related  to  the  zero-order 
r-matrix  by: 

(K,  K,- )  =  T(K)6(K  -  K,- )  (A7) 
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First-order 


The  first-order  7-matrices  are  proportional  to  the  2-D  Fourier  transform  of  the  surface, 
F(K),  and  can  be  written  as: 

7W(K,K,- )  =  Hn(K,K,)F(K  -  K,)  (A8) 

7/]\k,K,)=  Hi2(K,K,)F(K-K,)  (A9) 

where, 


F(K)  =  rf^J!/(R)e-*' « . 

(271:)^  •’ 

From  Moe  and  Jackson,  (1994)  Eqs.  (26)  and  (29) 
(i+r(K"))(i+r(K)) 


(A10) 


Hu{K",K)  =  - 


2/(3i(K") 


(All) 


i-1 


(i-hr(K"))(i+r(K)) 


VP  A  H 

Since  in  the  perturbation  expansion  (A1),  the  factor  jfej  is  displayed  explicitly,  the  above 
expression  is  actually  Eq.  (29)  of  Moe  and  Jackson  divided  by  a  factor  of  k^.  The 
following  expression,  equivalent  to  (All),  is  convenient  in  the  present  problem: 

1 


Hn(K",K)  =  - 


2iPl(K") 


{a(i+r(K"))(i+r(K))+i,(i-r(K"))(i-r{K))} 


a(K",K)  = 


--1 

VP  , 


K"  K  , 

- ;^+l - 


Z7(K",K)  =  Pi(K")Pi(K)(p-l). 


(A12) 

(A13) 

(A14) 
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From  Eq.  (24)  of  Moe  and  Jackson  (noting  the  difference  in  scalar  field  definitions,  and 
no  stratified  media,  ^23  =  0),  the  continuity  of  pressure  boundary  condition  for  first- 
order  fields  yields: 


(K")  =  -ij  d^K^i  (K)F(K"  -  K) 

x{iH^  1  (K",K)  +  (1  -  r(K))Pi  (K)  -  k(32  (K)(l  -h  r(K))} 

Another  expression  for  'F^|^(K")  follows  from  Eq.  (13)  and  (A9) 

'fW(K")  =  |j2f'F,(K)Hi2(K",K)F(K"-K). 

Equating  (A15)  and  (A16)  results  in 

H,2(K",K)  =  H„(K'',K)  +  (l-r(K))li,(K)  +  iKP2(K)(l  +  r(K)). 

Noting 

pPi  1+r 
kP,  i-r’ 

(A17)  simplifies  to: 

«,2(K",  K)  =  1  (K",K)  +  (1  -  p)|5,  (K)(l  -  r{K))} . 


(A15) 

(A16) 

(A17) 

(A18) 

(A19) 


Substituting  (A12)  into  (A19)  results  in  a  useful  expression  tor  Hi2(K",K),  given  as  Eq. 
(75). 


FIGURE  CAPTIONS 


FIG.  1.  Scattering  problem  geometry  for  rough  interface  separating  a  lossless  fluid  in 
medium  1  (z>/(R)),  from  a  lossy  fluid  in  medium  2  (z</(R)).  This  diagram  can  be 

viewed  as  a  slice  through  a  2-D  surface. 

FIG.  2.  Diagram  describing  geometry  of  variables. 

FIG.  3.  Illustration  of  effect  of  dispersion  on  peak  intensity.  Peak  of  g^{t),  normalized 
with  respect  to  as  a  function  of  grazing  angle. 

FIG.  4.  First-order  and  zero-order  field  strength  as  a  function  of  depth.  Field  strength 
given  in  dB  with  respect  to  incident  pressure.  Critical  angle  is  27.75°,  /  =  20kHz, 
a  =  0.1m,  7  =  3,  W2  =2xl0“^m,  p  =  2,  8  =  0.019,  q  =1500m/s,  v  =  1.13. 

FIG.  5.  Time-dependent  intensity  at  depth  0.2m  below  the  mean  surface.  Source 
frequency  =  20kHz,  /_j  =  100|ij.  Solid  line:  dispersion  terms  included;  dotted  line: 
dispersion  terms  set  to  zero,  g{t)  =  u{t).  Same  sediment  parameters  and  incident  angle 
as  Fig.  4. 

FIG.  6.  Zero-order  coherent  time  dependent  intensity  shown  with  dotted  line,  first- 
order  incoherent  intensity  shown  with  solid  line.  All  parameters  identical  to  Fig.  5. 
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